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Abstract 

Analysis of the electron density distribution in clusters composed of hydrogen flu- 
oride, water, and ammonia molecules, especially within the hydrogen-bond domains, 
reveals the existence of both a- and vr-binding between molecules. The cr-kind density 
distribution determines the mutual orientation of molecules. A vr-system may be de- 
localized conjugated, which provides additional stabilization of molecular clusters. In 
those clusters where the sequence of hydrogen bonds is not planar, a peculiar kind of 
TT-conjugation exists. HF2 and HsO^ ions are characterized by quasi-triple bonds be- 
tween the electronegative atoms. The most long-lived species stabilized by delocalized 
TT binding are rings and open or closed hoops composed of fused rings. It is conjugated 
vr system that determines cooperativity phenomenon. 

Keywords: molecular clusters, electron density distribution, hydrogen bonds, conjuga- 
tion, cooperativity 
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1 Introduction 

A fruitful idea abot the existence of a special kind of intermolecular bonds between such 
particles as water or ammonia was put forward already a century ago [Ij. Although T.S. 
Moore and T.F. Winmill [1] spoke only about a weak union formed by the molecules, and the 
concept of hydrogen bond developed later, the idea was already born. M.L. Huggins |2] and 
W.M. Latimer and W.H. Rodebush [3] are usually considered as the originators of the notion, 
but in its conventional variant it appeared and was accepted by the scientific community 
due to L. Pauling's [1]. Since that time, the scope of information about hydrogen bonds in 
the most diverse molecular systems and ensembles has kept increasing. As was mentioned 
in [5], if considering only the three- year period of 2006-2008, there appeared more than 
31000 papers in English, devoted to hydrogen bonding indexed in SciFinder. The amount of 
papers published in all languages for the last three quarters of a century would be impossible 
to estimate. It might even exceed a quarter of a million. 

It may be that no other problem of structure theory attracted so continuing attention 
of researchers. The reasons are quite clear. Hydrogen bonds dramatically change the prop- 
erties of molecular ensembles; and H-bonded ensembles are present nearly everywhere: in 
biological species, in the atmosphere, in many areas of technology. Plants and all living or- 
ganisms are composed of water, up to 90-99%. Tetrahedral arrangement of water molecules 
determined by the peculiarities of the electron density (charge) distribution was the first 
principal discovery in the field |6]. Later on, the scientific community realized that proteins 
exist in the form of a- helices and /3-sheets [2J, DNA molecules form double helices [8], and 
it is hydrogen bonds that stabilize these structures. Thus, all biological processes proceed 
in H-bonded environments and involve molecular structures with H-bonds. Insofar as the 
majority of technological processes proceed in aqueous solutions, and, which is more impor- 
tant, modern technology attempts to mimic the nature on the micro- and nanoscales, we 
can say that most of the technological schemes are also determined by peculiar properties of 
H-bonded substances. Finally, the formation of mist and clouds and the relevant processes 
in the atmosphere, both the troposphere and the stratosphere, are also determined by the 
ability of water molecules to aggregate either in the absence or in the presence of foreign 
particles. Thus, everything in and around us is related one way or another to the properties 
of H-bonded systems, but we still have no clear concept of H-bonding. 

2 New definition of hydrogen bond: progress and un- 
solved problem 

Having analyzed the milestone works in the field (refs. in [5]), members of the lUPAC Task 
Group on Categorizing Hydrogen Bonding and Other Intermolecular Interactions recom- 
mended in 2011 the following definition "The hydrogen bond is an attractive interaction 
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between a hydrogen atom from a molecule or a molecular fragment X-H in which X is more 
electronegative than H, and an atom or a group of atoms in the same or a different molecule, 
in which there is evidence of bond formation." X-H is the bond donor, while Y or Y-Z is 
the acceptor with an electron rich region, such as a lone pair of Y or a vr-bonded electron 
pair of Y-Z. 

The definition is accompanied by detailed comments and lists of criteria and character- 
istics of hydrogen bonds. The main features are summed up below [H]: 

(1) forces that determine the formation of a hydrogen bond have an electrostatic origin, 
arise from charge transfer between the donor and acceptor, and involve the dispersion 
component; 

(2) electron density distribution in the X-H. . . Y-Z domain is usually characterized by a 
bond path between H and Y with a (3, -1) critical point; 

(3) H-bond energy correlates well with the extent of charge transfer between the donor 
and the acceptor; and the strength of H. . . Y bond increases with an increase in the 
electronegativity of X; 

(4) X-H. . . Y angle is usually straight, and the closer it to 180°, the stronger the hydrogen 
bond and the shorter the H. . . Y distance; accordingly, hydrogen bonds are directional, 
which affects packing in the corresponding crystals; 

(5) X-H distance is larger when the fragment is involved in hydrogen bond, and the larger 
the distance, the stronger the H. . . Y bond; the H-bond formation is accompanied by 
the decrease in the X-H stretching frequency, the increase in the intensity of the 
corresponding infrared absorption line, and the appearance of new vibrations; 

(6) X-H. . . Y-Z hydrogen bond has characteristic NMR signatures including the proton 
deshielding, spin-spin coupling between X and Y, and nuclear Overhauser enhance- 
ments; 

(7) hydrogen bonds are involved in proton-transfer reactions (X-H. . . Y — X. . . H-Y) and 
may be considered as partially activated precursors of the reactions; 

(8) networks of hydrogen bonds show the phenomenon of cooperativity that is manifested 
in deviations of the properties from pair-wise additivity. 

The features are not included in the main definition body, because there are hydrogen 
bonds that possess not all, but only some properties. As examples, one can recall either 
blue-shifting X-H. . . vr bonds (X = O, N, C) [lOj or very weak F-H. . . Rg bonds (Kg = Ne, 
Ar) (which were classified as hydrogen bonds only due to the existence of an H. . . Rg bond 
path with a (3, -1) critical point) 
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Thus, all the characteristics of H-bonds, which were gradually discovered for the three 
quarters of a century, are accurately systematized now, but we still have no unique con- 
cept at out disposal. The concept is strongly needed, because until we clarify the nature 
of hydrogen bonds, we cannot progress in explaining the aforementioned geometrical and 
energetic properties of H-bonded systems and their experimental manifestations, especially 
the cooperativity phenomenon that underlies the functionality of all, primarily biological, 
species. Being unable to explain, we have no means for founded directional production of 
systems with desirable structures and properties, which is a serious obstacle on the way of 
developed nano technology. 

3 Two main kinds of hydrogen bonds 

H-bond energies depend on the nature of interacting molecules, their mutual arrangement, 
and coordination neighborhood and fall in a broad range, from 0.5 to 40 kcal/mol. The 
lowest energy is typical of Rg. . . HF complexes, while the highest is observed in H5O J (36 
kcal/mol) and HFg (40 kcal/mol) systems, which have symmetric configurations where a 
proton is located accurately at the midpoint between two oxygen or fluorine nuclei so that 
the structures are stabilized by two equivalent H. . . X bonds (X = F or O) rather than 
one covalent and one hydrogen bonds (X-H. . .X). Nuclear configurations of both molecular 
ions look as if these are unstable intermediates in the symmetric proton transfer reactions 
H3O+...H2O H2O. ..H3O+ or FH...F~ F"...HF. However, the configurations are 
stable when there is no external perturbation. When a perturbation appears, they do become 
unstable intermediates. A reason for that cannot be found if one bases the consideration on 
the above definition of hydrogen bond. 

Then, the following question arises: what is the kind of electron binding that determines 
the formation of hydrogen bonds? Is it reasonable to consider the interaction as no more than 
a perturbation of individual molecules, which is implicitly implied when one distinguishes 
different forces that contribute to the formation of hydrogen bonds? The approach is justified 
when one deals with typical intermolecular bonds of van-dcr-Waals kind, that is, when two 
approaching particles slightly affect the electron density distributions of each other, and the 
interaction can be estimated as a sum of the first- and second-order perturbation-theory 
corrections to the zero-order value that corresponds to non-interacting particles. The first- 
order contribution is the averaged electrostatic energy, while the second-order correction 
is the polarization (induction and dispersion) energy. Note that for estimating both first- 
and second-order energy corrections, it is sufficient to know the zero-order functions solely. 
However, being agglomerated in H-bonded clusters or ensembles, molecules no longer can be 
treated as independent individual particles, since they acquire novel properties; and hydrogen 
bonds formed (in most cases) have more in common with ordinary chemical bonds than with 
van-der-Waals bonds. Let us found and develop the idea. 

Note that molecules such as hydrogen fluoride, water, and ammonia, as well as their 
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so-to-say derivatives, namely molecules involving OH and NH groups, can form H-bonded 
networks. At the same time, there are examples of H-bonded molecular pairs that cannot 
be joined with each other. Therefore, it is reasonable to distinguish two kinds of hydrogen 
bonds: (i) separate bonds formed each by two particles and (ii) bond networks that unite a 
large number of molecules. The aforementioned weakest bonds are typically of the former 
kind, while the strongest bonds, in which FH, OH, or NH groups are involved, are of the 
latter kind. 

Bond networks appear when each involved particle can form no less than two bonds, 
one as a proton donor {d) and the other as an acceptor (a). This is possible when there 
is a proton, which is chemically bound to an electronegative atom, and a "free" valence 
electron density (uninvolved in interatomic bonds). The closer the proton and the "free" 
density domain in the particle, the stronger the correlation between hydrogen bonds. In 
other words, when hydrogen bonds are spatially close, i.e., separated with one intramolecular 
bond (X-H. . . X-H. . . X-H. . . ), a system of interacting bonds rather than just a spatial bond 
network originates. Below, we analyze electronic structures of the clusters composed of water, 
ammonia, and hydrogen fluoride molecules in order to clarify the character of the electron 
density distribution that determines the internally consistent state of molecular ensembles. 

Insofar as there is a large scope of data about cluster systems composed of water, am- 
monia, and hydrogen fluoride molecules, we believe it unreasonable and impossible to quote 
all of works and, at the same time, disrespectful to quote only some of them and neglect the 
residual. Therefore, the discussion below contains no reference to independent studies. A 
representative list of references can be found in the aforementioned technical report of the 
lUPAC Task Group [5j. 

4 Electron density distribution: CO LCMO approach 

Undoubtedly, molecular orbital is a conditional mathematical construction. Nevertheless, 
electrons are physical objects, and spatial distributions of their charges can be described by 
certain functions. Until we clarify the physical nature of the objects, we can construct only 
approximate functions that represent (more or less accurately) the mean exterior of the real 
charge distribution. Then, at present moment, the most reasonable approach is to consider 
natural orbitals, in the basis of which the total electron density of a molecular system p{r) 
is written as a sum of the weighted densities that correspond to individual 0i(r) orbitals: 

p 

i=l 

When Hi weights (or coefficients), which are referred to as effective occupancies of the or- 
bitals, are close to one or two (i.e., negligibly differ from these integers), one can believe 
(just believe, since this cannot be proved) that 0i(r) functions, which were found as a re- 
sult of nonempirical calculations, represent reasonable approximations of the actual charge 
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distributions of individual electrons. 

In the case of water, ammonia, and hydrogen fluoride molecules (in the ground elec- 
tronic states), as well as clusters composed of these molecules, rii numbers in the above 
equation meet this requirement; and natural orbitals constructed in the second-order of the 
Moeller-Plesset perturbation theory (MP2), as well as with the use of fully or partially vari- 
ation schemes, such as single- or multi-reference configuration interaction method or coupled 
clusters approach, are very close to the solutions obtained at the Hartree-Fock (HF) level. 
This means that if we were interested only in the orbital picture of the systems, we could 
restrict the consideration by the HF level. However, when estimating relative stabilities 
of diverse cluster configurations or dissociation energies, we cannot neglect the correlation 
energy contributions. 

Note that it is MO LCAO scheme (where molecular orbitals are approximated by lin- 
ear combinations of atomic orbitals) that enables one to propose the most straightforward 
physical interpretation of the quantum chemical description of the electronic subsystems 
of individual molecules. It is founded when a single-determinant description of the elec- 
tronic system provides reasonable results. Let us use a similar approach for interpreting the 
electronic wavefunctions of molecular clusters implying that cluster orbitals (CO) can be ap- 
proximated by linear combinations of molecular orbitals (LCMO). Such CO LCMO scheme 
is no less founded than conventional MO LCAO one, since relative changes in the electron 
density distribution of individual molecules that accompany the formation of intermolecular 
bonds are definitely smaller than those typical of intramolecular chemical bonds. 

Results considered below are obtained in quantum chemical calculations carried out at 
the MP2 level with a 6-3 IG basis set augmented with diffuse and polarization functions on 
all nuclei (6-31++G(d,p)). This level of theory was repeatedly proved to be adequate for 
describing the systems of interest. All the configurations correspond to the energy minima, 
which is confirmed by the normal-coordinate analysis. Dissociation energies are estimated 
taking into account counterpoise corrections for the basis set superposition error, and the 
mean H-bond energy is found by dividing the total dissociation energy of a cluster (X„ 
—7- nX, where X = H2O, HF, or NH3) by the number of H-bonds. All calculations were 
carried out with the use of GAMESS [12] and Firefly [I3] program packages. 

4.1 HF, H2O, and NH3 molecules: MO basis 

Orbitals of individual water, ammonia, and hydrogen fluoride molecules, which serve as 
the basis at the stage of constructing cluster orbitals, are well known and look as quite 
simple combinations of atomic orbitals. Let us denote orbitals of individual molecules as 
(p^, and 0° in the case of hydrogen fluoride, water, and ammonia respectively, while 
those of clusters as ip{ , ipf, and 93^ u. Being strongly localized Is core orbitals of fluorine, 
oxygen, and nitrogen atoms respectively, 0^, and 0" molecular orbitals are not of interest 

^ Sometimes, superscripts are omitted when no confusion may arise (e.g., in figures). 
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from the point of view of the formation of intermolecular bonds. Below, we briefly remind 
how the bonding orbitals of individual molecules look like omitting in the mathematical 
expressions those atomic orbitals, which contribute with negligible or much smaller weights, 
and dropping coefficients that depend on internuclear distances and particular functional 
form of the orbitals. 

In an HF molecule (Fig. chemical bonding is determined by two orbitals, one of 

which (61) is a linear coniination of 2.(F) and ^(Ha aton.. taction, while the other 

{4>l) is a combination of 2pz(F) and ls(H) functions {Oz axis is oriented along the F-H 
line). The residual four valence electrons are described by two degenerate orbitals, namely, 
(f){ = 2px{F) and (pi = 2py{F), and represent the electron density that can be referred to as 
"free", i.e., uninvolved in interatomic binding. 



Figure 1: Valence molecular orbitals of individual (a) HF, (b) H2O, and (c) NH3 molecules. 

In an H2O molecule (Fig. Wp), the lowest valence orbital ((/)^) is a linear combination of 
2s(0), Is(H^), and Is(H^) atomic functions (where superscript enumerates protons). Two 
next orbitals are bonding in and H^. . . domains respectively: 0^ = 2^^.(0) 

+ (Is(H^) - Is(H^)) {xOz is the nuclear plane of the molecule, and Oz axis is the symmetry 
axis of the molecule) and (p^ = 2pz{0) + (Is(H^) + Is(H^)). Finally, a lone pair, which 
represents the free density of the molecule, is described by (p^ = 2py{0) orbital. 

^Figures [U [SHHl [3 and IMTTI show the surfaces of constant p — 0.0015 01 ip ^ 0.0015 value, i.e., p = 
2.25-10^^ au. The darker regions correspond to the positive values oi p ox ip function, while the lighter 
regions, to the negative values, so that their boundaries are the zero electron density surfaces. Orbitals in 
Figs. [1] and [3] are shown as transparent so that one can clearly see positions of nuclei. 

■^We use conventional notations for the atomic orbitals that involve shell numbers. 



^1 



4>3 



<1>4 



^5 
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In an NH3 molecule (Fig. Wp), the lowest bonding orbital is (j)2 — 2,s(N) + (Is(H^) + 
Is(H^) + ls(H'^)). Two next degenerate orbitals (^3 and (pD provide the electron binding 
in the N-H\ N-H^, and N-H^ domains: 0^ = 2p^(N) + (1s(H2) + ls{B.^) - ls{}i^)) and 
04 = 2py{N) + (Is(H^) - Is(H^)) {Oz axis coincides with the triad axis of the structure, 
and Ox axis coincides with the projection of N-H^ bond on the orthogonal plane). The 
residual valence orbital is also bonding: 0^ = 2pz(N) - (Is(H^) + Is(H^) + Is(H^)). Thus, 
by contrast to water and hydrogen fluoride molecules, there is no electron uninvolved in 
intramolecular bonds in ammonia; only about a half of the electron density described by the 
latter orbital and localized outside the NH3 pyramid can conditionally be referred to as free. 

Thus, if one supposes that the main contribution to the formation of hydrogen bonds 
is provided by the redistribution of the "free" electron density of the molecules, the most 
spatially flexible systems of H-bonds with the highest bond energies should be observed in 
the case of HF molecules and the smallest flexibility and the weakest bonds, in the case of 
NH3. Let us check the supposition by considering some typical small (HF)„, (H20)„, and 
(NH3)„ clusters. 

However, before proceeding to particular examples, it is worth noting that intramolecular 
bonds form due to the overlap of the frontier orbitals of interacting particles (atoms or func- 
tional groups). These are both highest occupied and lowest unoccupied molecular orbitals 
or singly occupied molecular orbitals of both reagents; and in the resulting system (in the 
ground electronic state), bonding orbitals are usually occupied, while antibonding ones are 
unoccupied. When one deals with clusters composed of closed-shell molecules, the situation 
is different: the redistribution of the electron density is due to the overlap of the occupied 
orbitals of the particles, so that the resulting functions, which describe the electron density 
distribution in the cluster, are of both kinds, namely bonding and antibonding. Hence, the 
total binding effect is much smaller compared to typical interatomic bonds and depends on 
the balance between bonding and antibonding contributions in each intermolecular region. 

4.2 (HF)„ clusters 

HF molecule typically forms two bonds with the neighbors {da coordination), though there 
are situations when either hydrogen or fluorine atom or both of them are involved in two 
bonds with the neighboring molecules. Stable structural units are molecular rings (Fig. |2^) 
where molecules are arranged in such a way that H. . . F-H angles are close to 90°, while 
dihedral F-H. . . F-H angles are arbitrary. When no less than six molecules are joined in 
a ring, it looks like a folded sheet or strip composed of rectangular segments (Fig. [2)d). 
Openwork species, in which three-membered molecular rings have common HF edges, are 
also possible (Fig. [Jt). Larger rings that involve three-coordinated molecules can be joined in 
nonplanar openwork species. Such building blocks enable one to construct a two-dimensional 
network, but any attempt to build a continuous three-dimensional bonded system fails. 
Whenever relatively small ring-like fragments can separate, they do. Let us try to understand 
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the peculiarity by considering bond energies typical of different coordination neighborhoods. 




Figure 2: Typical structures of (HF)„ clusters: (a) plain rings; (b) extended folded rings; 
and (c) openwork species composed of trimolecular rings. 

When the number of molecules in (HF)„ rings (Fig. [2^) increases from 3 to 5, the mean 
H-bond energy changes from 5.2 to 8.6 kcal/mol. In extended folded rings similar to that 
shown in Fig. |2]3, a mean bond energy is about 8.4 kcal/mol. Two parallel molecular 
rings have no noticeable attraction to each other. Species based on trimolecular rings with 
common HF edges (Fig. ^) are less stable than plain rings: a mean H-bond energy is about 
4.2 kcal/mol. With an increase in the number of molecules, which form more than two 
hydrogen bonds, the mean bond energy decreases to about 4 kcal/mol. A mean bond energy 
per molecule is higher, about 5 kcal/mol, but still substantially lower than that in plain 
molecular rings. Bonds in chains are characterized by intermediate energies, e.g., about 5.9 
and 6.4 kcal/mol in four- and five-membered chains. 

Thus, knowing the bond energies in diverse structural units, one can predict most prob- 
able kinds of H-bonded systems. However, the information cannot clarify the underlying 
reasons of stability. For example, the difference between the bond energies in molecular 
rings and chains (2-2.5 kcal/mol) seems too large for the systems composed of molecules 
that have no more than two neighbors. At the same time, when the number of H-bonds 
formed by a molecule increases to three or four, their energies decrease despite nearly the 
same change in the atomic charges, which a re conventionally thought of as correlating with 
the bond strengths, hence, the bond energies depend on the character of the electron density 
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distribution rather than its integral value in the corresponding spatial regions. Peculiarities 
of the distribution can be judged from the character of cluster orbitals. 

Let us start with (HF)3, which is the smallest ring-like structure. Even a prompt glance 
distinguishes two kinds of cluster orbitals (Fig. [3]). The lowest Lpl^ orbitals with m = 4, 
5, and 6 (three core functions are again omitted) are obviously linear combinations of 02 
molecular orbitals and have definitely a cr-kind. To the approximations of all the residual 
cluster orbitals, p-functions of fluorine atoms noticeably contribute. Six of them, namely, 
(fi{^ with m = 7, 8, 9, 11, 14, and 15, are also of a-kind and describe either electronic 
binding or antibinding in the H-bond domains. Let us denote the plane of fluorine nuclei 
as xOz; then, these cluster orbitals can be considered as linear combinations of and 
molecular orbitals. Note that the residual cluster orbitals with m = 10, 12, and 13 have 
definite 7r-character and are approximated by linear combinations of (f)( molecular orbitals. 
These cluster orbitals have much in common with typical 7r-system functions of a conjugated 
hydrocarbon. The sole peculiarity is as follows. In hydrocarbon rings, maxima of the vr 
functions are at carbon nuclei, while minima and zeros are at the midpoints of C-C bonds. 
In (HP) 3 cluster, the large- value parts of orbitals correspond to H-P regions with maxima at 
P nuclei and minima in H. . . P domains close to bridge proton positions. Thus, the orbitals 
describe a delocalized ir-binding of fluorine atoms; and H atoms here act as auxiliary bridge 
sites. 

All a and tt cluster orbitals remain bonding in the regions of intramolecular H-P bonds, 
but can be either bonding or antibonding in the hydrogen-bond domains. Por example, of the 
three cluster orbitals that can be constructed (taking into account the general C3 symmetry 
of the cluster) from three (jy^^^^ molecular orbitals {i superscript enumerates molecules), only 
one {ip{) is bonding, while two residual have prevailing antibonding character: 

= 0f ) + 4''' + 4''' 
4 = 20j(^) - 4^'^ - 

Insofar as a cluster orbitals are linear combinations of either (p^ = 2s(P) + ls(H) or 
03 = 2pz{F) + Is(II) bonding orbitals and 0{ = 2px{F) lone pair orbitals, whose axes 
are mutually orthogonal, HP molecules tend to be arranged in clusters in such a way that 
H. . . P-H angles are close to 90°. Planar configuration of (HP)3 ring also seems natural, 
since the additional tc binding is provided by the overlap of (p( = 2py{F) molecular orbitals 
whose axes are normal to the plane of cr-binding. 

The same character is typical of the four-membered (HP)4 ring and openwork (HP)„ 
clusters, in which trimolecular rings are joined via common HP edges (see Pig. [2]). On the 
whole, changes in the orbital picture that take place with an increase in the ring size or upon 
fusion of the rings are similar to those observed in the case of hydrocarbon rings when the 
number of C nuclei in the ring skeleton increases or when rings of the same size are joined 
via edges. An interesting illustration is given by Pig. H] where five occupied tt orbitals of 
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Figure 3: Valence orbitals of (HF)3 cluster. 

(HP) 5 openwork structure are compared to five tt orbitals of naphthalene. In the case of 
fused rings, the electron density, which is "contributed" by a common edge HP molecule 
to the formation of hydrogen bonds, is nearly twice as small as the density, which provides 
binding of molecules in plain rings. As a result, the energies of such bonds are noticeably 
lower. 

The similarity of the orbital pictures of small (HP)„ rings and aromatic hydrocarbons 
does not apply to chain-like structures. In planar molecular (HP)4 or (HP)5 chains, there 
are both a and tt cluster orbitals, but each tt orbital describes chiefiy the electron binding 
of some two neighboring molecules. Thus, by contrast to rings, tt orbitals in the chains are 
noticeably locahzed. This is clearly refiected in the internuclear distances. Por example, 
in (HP)4 chain, the increase in H-P distances in the central molecules (compared to the 
individual molecule) is about 0.02 A while that in the terminal molecules is no more than 
0.01 A. Internuclear distances are not equal and fall in a range from 1.636 to 1.714 A. At 
the same time, in the ring of the same molecular size, the increase in the intramolecular 
H-P distances is 0.035 A, and the intermolecular H. . . P distances are all the same, about 
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1.534 A. 

The more localized n bonding orbitals, the weaker the bonds; and the difference in 
energy is substantial. When a zigzag tetra- or pentamolecular chain closes in a ring, the 
absolute energy decrease is 17-20 kcal/mol, which is nearly thrice as large as an individual 
H-bond energy in the ring. This means that such closing is accompanied by the additional 
stabilization of the already existing hydrogen bonds, which is reached due to the appearance 
of a delocalized conjugated vr system. 

To sum up, the combination of a and vr binding explains the relative stability of small 
planar (HF)„ clusters. Planarity is the characteristic feature of conjugated hydrocarbons. 
However, as was mentioned above, when there is no less than six molecules in an (HF)„ ring, 
it is no longer planar (even at n=5, there is already a small deviation from planarity). It is 
folded and consists of rectangular segments. This peculiarity agrees with the trend of HF 
molecules to form sequences of mutually orthogonal hydrogen bonds, which is impossible in 
the case of planar structures. The orthogonality is determined by a binding. Hence, we can 
state that this kind of electronic binding is present in all the clusters and determines the 
spatial arrangement of molecules. Then, what is the role of n orbitals and do they exist in 
nonplanar systems? It turns out that they do, but in some new variant, which is unusual 
for any person who has some acquaintance with organic chemistry. 

As an example, let us consider tt orbitals of a ring-like (HF)6 cluster (Fig. The lowest 
bonding orbital of this kind [ipa) has much in common with usual tt orbitals of conjugated 
hydrocarbons with a sole difference: with an increase in the number of molecules, minima 
of the function (more accurately, (3, -1) points) become closer to the position of the bridge 
proton. (A similar trend is observed in the case of a cluster orbitals.) All the residual vr 
cluster orbitals dramatically differ from those of conjugated hydrocarbons and planar (HF)„ 
clusters. Due to the planar configuration of all the latter systems, there were definite charge- 
alternation surfaces in all of them located in the H-bond domains. 

In (HF)6 cluster, constant-sign parts of the orbitals look as closed rings, and the rings 
are connected as links of a chain {(ff, in Fig. [S^). The highest orbital resembles a chain rivet 
{ipc in Fig. Et), as if in a six-link ring of one sign, a rod of the other sign is inserted, and 
both its ends are riveted to triangles. Thus, the antibonding regions are located either in 
the vicinity of bridge protons or in the central part of the cluster. And electrons of all six 
TT orbitals, each of which is delocalized over the whole hexamolecular ring, contribute to the 
intermolecular binding. 

When the number of molecules is increased above eight, and pronounced folded structures 
appear, the orbital picture becomes more complex (Fig. [5]d). Even the lowest tt orbital looks 
already as two interlaced rings (of different signs) and resembles a typical tt bonding orbital 
only in a two thirds of the ring {(pd)- In the next orbitals, which are more and more complex 
interlacements of multifold bent rings, the sign-alternation regions already coincide with the 
H-bond domains, so that the electronic binding in certain H. . . F regions is weaker and, at 
the same time, there is no general pronounced binding of all the molecules. 
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a b 




Figure 5: Typical valence 7r-orbitals of nonplanar (HF)„ clusters: (a) n = 6 and (b) ?7, = 9. 

In a general case, if we look at a folded structure from the side, we can see a zigzag; and 
the more strongly internally bound subsystems correspond to separate periods of the zigzag. 
For example, in (HF)9, tetra- and pentamolecular subsystems {ipe and (ff in Fig. Eb) can be 
distinguished. Quantitatively, this separation of quasi-aromatic segments in extended folded 
structures is reflected by the mean H-bond energy of the cluster, which decreases from 9.2 
kcal/mol in (HF)6 to 8.3 kcal/mol in (HF)9. 

The discovered peculiarities of electronic binding between HF molecules show that hydro- 
gen fluoride should form structures composed of lamellar domains, the most stable building 
units of which should be planar or folded (depending on the molecular size) rings. In the 
rings, molecules are joined via electronic bonds, in which a and tt contributions are clearly 
distinguishable, and vr binding system is delocalized conjugated. With an increase in the 
molecular size of the ring, which is accompanied by the formation of nonplanar structures, 
the H-bond energies increase (to a certain limit) due to the stronger vr binding. When the 
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ring becomes too large, and no effective conjugated system of bonds can be formed, strongly 
bound subsystems are separated. 

Note tliat two parallel zigzag chains located at a distance of about 2.5 A (equal to F- 
H. . . F distance) should form a similar vr-bonded agglomerate. However, in the absence of 
terminal closing molecules, which keep the chains at a constant distance, such agglomerate 
should exist only at a low temperature. With an increase in temperature, it should probably 
transform in several rings arranged successively. At the same time, antiparallel zigzag chains 
can easily be interconnected via properly reoriented molecules, so that either one extended 
ring or a similar succession of rings can appear. 

These peculiarities explain why hydrogen fluoride is a mobile and volatile liquid. 

Thus, in any case, we deal with a coherent binding of all the molecules or some their large 
subgroups rather than simply binding of molecular pairs. Let us now see whether a similar 
conjugation exists in clusters composed of water and ammonia molecules, the ensembles of 
which cannot have planar configurations, but X-H. . . X-H. . . bond sequencies (X = O or N) 
can form planar rings in some cases. Note that both H2O and NH3 molecules can form up 
to four "normal" hydrogen bonds with the neighbors, which makes them more promising as 
building blocks of three-dimensional networks. Let us start with ammonia clusters, which 
should be in some sense opposite to hydrogen fluoride ones, since the "free" electron density 
of ammonia molecules is the lowest (compared to HF and H2O). 

4.3 (NHs)^ clusters 

Ammonia molecule can normally form three hydrogen bonds as a proton donor and one as 
an acceptor. When its coordination neighborhood is of da kind (Fig. [H^), i.e., molecules 
form a ring-like structure, the mean bond energy is nearly by half smaller compared to that 
in (HF)„ clusters. For example, when the number of molecules in a ring increases from 
three to six, the energy increases from 3.6 to 4.5 kcal/mol. When small molecular rings are 
joined in branched or cage-like structures via molecules whose nitrogen atoms are involved 
in two hydrogen bonds with the neighbors (Fig. Eb), the mean bond energy is lower: about 
3.1 kcal/mol. In such systems, similarly to (HF)„ clusters, the bond energy per molecule is 
higher (about 3.8-4.0 kcal/mol), but still lower than that in plain rings. 

The formation of hydrogen bonds between ammonia molecules is accompanied by such 
electron density redistribution that practically does not touch those hydrogen atoms, which 
are not involved in the bonds formed. The extent of integral charge redistribution is com- 
parable to that in (HF)„ clusters, but the bond energies in the systems differ nearly by half. 
This indicates again that integral electron density characteristics can by no means be unique 
indicators of the H-bond strength. 

Orbital picture of (NH3)3 is similar to that obtained in the case of hydrogen fluoride 
trimer. Twelve valence cluster orbitals (v?^) can be divided in four threes (m = 4, 5, 6; 7, 11, 
12; 8, 9, 10; 13, 14, 15), each of which is represented by linear combinations of the molecular 
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Figure 6: Typical structures of (NHa)^ clusters: (a) plain rings and (b) joined or fused tri- 
and tetramolecular rings. 

orbitals of the same kind, namely, 03; 4^2^ respectively. Here, 0^ molecular orbitals 

with A; = 2, 3, or 5 provide combinations of cr kind, whereas overlapped 02 orbitals give tt 
functions, one of which {ip'^) is a typical bonding orbital, while two residual (999 and (/?"q) 
are degenerate and describe electron bonding of some two molecules and antibonding in two 
residual pairs. 

Thus, (7- binding of ammonia molecules is provided by the "free" electron density of the 
molecules (0f functions) and the electrons described by one of two degenerate molecular or- 
bitals (^3). Electrons that were described initially by 02 orbitals (orthogonal to functions) 
represent the conjugated tt system of the molecular ring. Orientation of 05 orbital along the 
triad axis of NH3 pyramid determines the preferable linear configuration of N. . . H-N frag- 
ments, i.e., the linearity of hydrogen bonds; whereas the original degeneracy of 03 and 02 
orbitals of individual molecules leads to the arbitrary orientation of N(^)H3 pyramid with 
respect to N^^^-H hne in (H2)N(^)-H. . . N(^)(H3) hydrogen-bonded pair and an arbitrary rota- 
tion angle of N^^^Hs pyramid about its triad axis. Accordingly, possible mutual orientations 
of NH3 pyramids in (NH3)„ clusters can be very different. The relative closeness of HNH 
valence angle in the individual ammonia molecule (ca. 113°) to the planar angle in an equi- 
lateral triangle provides the formation of such (NH3)3 structure where (N-H)3 bond ring is 
planar and have much in common with the (HF)3 ring. The difference is the substantially 
lower electron density in the H-bond domains, which is determined by the fact that it is 
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electrons, which provide the existence of intramolecular N-H bonds, rather than lone pair 
electrons that contribute both to a and tt binding. As a result, bond energies are much 
smaller: 3.5 kcal/mol in (NH3)3 versus 5.3 kcal/mol in (HF)3. 

In a tetramolecular ring, the orbital picture is principally similar; though, due to a small 
deviation of (N-H) 4 ring from planarity, nonbonding and antibonding tt orbitals are already 
slightly similar to interlaced rings of different sign. However, the electron density distribution 
within such rings is substantially nonuniform, with constrictions in the H-bond domains close 
to proton positions. 



a b 




Figure 7: Typical valence tt orbitals of (NHs)^ clusters: (a) n = 6 and (b) n = 5. 

In the case of a hexamolecular (NH3)6 ring, which has a bath-like configuration, vr-orbitals 
look like constant-sign rings interlaced in a complex way (Fig. [7^); and the rings are broken 
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except for the sole orbital, ipa- The latter statement is conditional, since in a const ant- density 
map with a boundary value smaller than 10~^ au, the rings will be seen as closed, but such 
low electron density provides actually no binding in the corresponding spatial regions. As 
can be seen, delocalized vr bonds exist within open three- or five-membered rings. 

A similar picture is observed in (NH3)5 cluster when it has a configuration of two tri- 
molecular rings joined via a common molecule (Fig. Wp)- The latter molecule has a ddaa 
coordination, so that despite its low "free" density it acts as an acceptor of two protons. Even 
the lowest vr orbital of the cluster [ipd in Fig. [7)d) describes a noticeable electron binding of 
only four molecules, i.e., except for the tetracoordinated central one. The residual vr cluster 
orbitals determine binding within either dimeric fragments or a dimer and a chain-like trimer 
{iPe and iff 'va. Fig. Wp)- In all the cases, the electron density within the spatial region of the 
central molecule and those its hydrogen bonds, in which it acts as an acceptor of protons 
(i.e., an electron density donor), is the lowest. 

Thus, the very low electron density, which can be contributed by an ammonia molecule to 
the formation of tt system can provide additional stabilization of only small rings (composed 
of three to five molecules) or small chain-like fragments within larger clusters. In the regions 
of tetracoordinated molecules, there is no continuous delocalized vr system that could stabilize 
the joint. Therefore, it is highly probable that the structure reorganization in such regions 
can be initiated by a very small external perturbation (or energy fiuctuation). The result will 
be the breakage of one or two H-bonds formed by the central molecule, and the appearance 
of bent chain-like fragments only partially stabilized by tt binding. 

Ring-like species with quasi-aromatic vr systems, which were found to be most stable in 
the case of (HF)„ clusters, are very small (three- to four-molecular) and statistically rare 
objects in ammonia molecular ensembles. This peculiarity explains, particularly, the much 
lower boiling point of liquid ammonia compared to hydrogen fiuoride. 

Now, let us turn to water clusters whose molecules have free electron density, which 
is noticeably larger than that of ammonia molecules, but smaller than that of hydrogen 
fluoride. 

4.4 (H20)„ clusters 

The mean bond energy in (H20)„ clusters is known to increase with an increase in n up to a 
certain limit, which is reached at a number of molecules about twenty. Here, people typically 
speak about the bond energies per molecule. However, the interaction energy per H-bond 
taking into account the neighborhood of molecules is more informative. Let us distinguish 
the following coordination kinds: da, dda, daa, and ddaa. If all molecules are involved in 
trimolecular rings, then, practically irrespectively of their particular local neighborhoods {da 
in a ring-like trimer or da and ddaa in a pentamer where two trimolecular rings are joined 
via a common molecule. Fig. [8^), the mean H-bond energy is 5.0 kcal/mol. In (H20)„ rings 
of the larger size (Fig. [8)d) where H-bond angles are closer to 180° and all molecules have 



18 



da coordination, the interaction energy per bond is 6.5, 6.8, and 6.9 kcal/mol at n = 4, 5, 
and 6 respectively. In (H20)6 prism, (H20)8 cube, and (H20)i2 cage, which are composed of 
fused three-, four-, five-, and six-membered molecular rings where tricoordinated molecules 
prevail (Fig. |8^, c), a mean H-bond energy depends on the steric factor, being 4.7, 5.7, and 
6.0 kcal/mol, respectively. Thus, with an increase in the number of hydrogen bonds formed 
by the same molecule, the mean bond energy slightly decreases, while with an increase in 
the size of the ring, it increases, so that the H-bond energy per molecule increases with 
an increase in the molecular size of the cluster, being 6.9, 8.5, and 8.6 kcal/mol in the 
aforementioned (H20)„ clusters with tt, = 6, 8, and 12. 




Figure 8: Typical structures of (H20)„ clusters: (a) composed of or involving trimolecular 
rings; (b) plain rings with nearly straight hydrogen bonds; and (c) cages composed of fused 
tetra-, penta-, and hexamolecular rings. 

Thus, water molecules seem to be most promising for the formation of three-dimensional 
bond networks. Energies of individual bonds in plain ring-like (H20)„ clusters are higher 
than those in (NHa)^ clusters and lower than those in (HF)„. However, it is water clusters 
solely where the mean bond energy remains sufficiently high when the coordination number 
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of the molecule increases from two to four; and the H-bond energy per molecule even increases 
(up to a certain limit). 

By contrast to both ammonia and hydrogen fluoride, water molecule has no degenerate 
orbitals. Recall that there are two degenerate orbitals in HF that describe the electron 
density of lone pairs; and there are two degenerate orbitals in NH3 molecule that describe 
N-H binding. When molecules of both latter kinds form clusters, electrons of one of two 
degenerate orbitals contribute to a binding, while those of the other, to the formation of 
TT-system. Orbitals of H2O molecule are all nonequivalent in energy and electron density 
distribution. Therefore, the orbital picture should strongly depend on the arrangement of 
molecules, namely, on the relative orientations of their nuclear planes. 

Nevertheless, based on the orbital picture of individual water molecule and the discovered 
peculiarities of the intermolecular bonds in (HF)„ and (NHs)^ clusters, we can expect that 
a binding in (H20)„ clusters should chiefly be provided by the overlap of intramolecular 
bonding orbitals, namely, c/)^ = 2p^(0) + (Is(Hi) - ls{}i^)) and 0^ = 2p^(0) + (Is(Hi) + 
Is(H^)), while lone pair electrons should provide n binding. It is quite possible also that lone 
pair electrons can be involved in the formation of a bonds, while the overlap of intramolecular 
bonding orbitals can additionally contribute to the formation of tt system. In the molecule 
that acts as a proton donor, electrons described by and ^4 determine a binding. In the 
proton-acceptor molecule, diverse MOs can contribute to a binding depending on the mutual 
orientation of the nuclear planes of the molecules. Such ambiguity makes it more difficult 
to predict optimal mutual orientations of water molecules in clusters and the preferable 
coordination of molecules; but, at the same time, it determines the larger flexibility of the 
electron density distribution between a and vr systems. 

Let us start (as before) with the smallest ring-like cluster, namely trimer. Recall that both 
in (HF)3 and (NH3)3, we could definitely distinguish a and tt cluster orbitals, and the latter 
orbitals were very similar to bonding and antibonding orbitals of a conjugated hydrocarbon. 
What do we see in the case of water trimer? Figure [9^ shows three most interesting valence 
orbitals of the cluster. Their approximations involve (p^ and (p^ molecular orbitals and they 
look like twisted tt orbitals. They are twisted in such a way that within each O. . . H domain 
they have definitely tt kind. Alternation-sign surfaces of the orbitals always pass in close 
vicinity to bridge protons. 

The orbital picture of the ring-like water tetramer also represents an intermediate situa- 
tion between those typical of planar and non-planar systems (Fig. [Hb). There is an orbital, 
which resembles a nonbonding tt orbital of a planar four-membered ring {(fia)- Another one 
{ipb) is a TT-bonding orbital. Two residual orbitals are similar to those observed in non-planar 
clusters of hydrogen fluoride and ammonia, namely, look like interlaced constant-sign rings 
(Fig. [Hb shows one of them, ipc), where the sign- alternation surfaces pass in close vicinity to 
bridge-proton positions. It is worth noting that, despite the visible resemblance, the orbitals 
have a definite peculiarity. Zero-density surfaces are very complex, and in many cases, they 
divide two bonding domains, namely, intramolecular 0-H and the neighboring intermolecu- 



20 




Figure 9: Some valence vr orbitals of ring-like (H20)„ clusters: (a) n = 3 and (b) n = A. 

lar H. . . O domains, so that the generally 7r-kind orbital that spans the whole cluster has a 
mixed vr/cr character in the H-bond domains. 

In ring-like pentamer and hexamer, the situation is principally similar. There are five 
and six cluster orbitals of generally tt kind respectively. Character of the electron density 
distribution described by any of the orbitals can be named intermediate between those 
typical of hydrogen fiuoride and ammonia clusters. This means that there are bonding 
and antibonding domains, but the constant-sign rings always look like bended or twisted 
doughnuts with tight constrictions. In fact, this observation strongly correlates with the 
"amount" of free electron density that can be localized in the H-bond domains. 

This result explains also why the mean bond energies in ring-like water clusters are 
intermediate between those in hydrogen fiuoride and ammonia clusters of the same size and 
why the difference in bond energies between hydrogen fiuoride and water clusters increases 
with an increase in the number of molecules involved in the ring. It is about 0.1, 1.0, 1.6, 
and 2.4 kcal/mol in the case of trimers, tetramers, pentamers, and hexamers, respectively. 
Thus, the most stable are rings composed of HF molecules. 

However, the situation dramatically changes when one goes to cage-like clusters. HF 
molecules do not agglomerate to form stable three-dimensional cages. By contrast, water 
molecules form three-dimensional networks. Mean H-bond energies in such systems are 
close to or even higher than 6 kcal/mol, and the bond energies per molecule are close to 
8.5 kcal/mol., i.e., of the same order as in the most stable folded (HF)„ clusters. The 
reason becomes clear if we consider Fig. [TOl where orbitals of (H20)5 cluster, which involves 
a tetracoordinated molecule, and those of (H20)8 and (H20)i2 cages composed of fused 
molecular rings, are shown. One of the orbitals of (H20)5 {^pa in Fig- ITUb) looks like two 
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concatenated H letters, one leg of each being normal to the residual T-shaped segment. Such 
orbital provides 7r-electron binding in the domains of all six H-bonds. Two other orbitals 
(nearly degenerate) look like concatenated twisted Z-shaped patterns (Fig. [TOk shows one 
of them, ifib) and provide n binding within one structural triangle and between two residual 
molecules. Similarly to the above orbitals of ring-like clusters, these orbitals have generally 
TT character and a mixed n/a character in the H-bond domains. Such flexibility of orbitals is 
determined by the contributions of ^3 , 04 , and (pf functions with different weights. Thus, 
a relatively small cluster that involves a tetracoordinated joint molecule is stabilized by vr 
binding between all molecules. Recall that, in the case of ammonia molecules, such binding 
was impossible. 

b 







Figure 10: Some valence tt orbitals of cage-like (H20)n clusters composed of fused or joined 
molecular rings and involving tri- and tetracoordinated water molecules: (a) n = 5, (b) n = 
8, and (c) n = 12. 



Orbitals that span the whole or large part of the structure can be found also in (H20)8 
and (H20)i2 cages. Figure fTOl (b. c) shows four such orbitals. In octamer, there are orbitals 
of diverse shapes, and here, the principally new situation is observed. First of all, there is 
an orbital {(pc), which looks like a sandwich and describes conjugated vr-binding systems of 
two four-membered rings that form two opposite faces of the cube (let us conditionally refer 
to them as the bottom and top). Another orbital of the same kind {ipd) describes binding 
within four tetramolecular rings that form side faces of the cube. There are also orbitals 
that describe tt binding in the domains of dimeric edges of the cube. Note that in all the 
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cases, the orbitals have bonding or antibonding a character in the directions normal to tt 
bonds. This is the pronounced manifestation of the aforementioned flexibihty of the electron 
density distribution in water clusters. 

In the case of (H20)i2, let us consider also only two orbitals (</?e and (/?/). At first glance, 
they are very similar to those of plain rings. In fact, they are. In one of them, interlaced 
constant-sign domains are composed of six- and five-membered rings connected via a typical 
TT bonding segment of a tetramolecular ring. In the other orbital, two concatenated constant- 
sign patterns are based on a six-membered ring fused with five- and four-membered rings on 
opposite edges. As a result, there appears a conjugated tt system that spans four structural 
rings, which in turn form a closed or open spatial hoop. Some other cluster orbitals look like 
interlaced or concatenated variously bent segments, and the corresponding electrons provide 
binding within four- to eight-molecular subsystems of the cluster. 

These results explain why cage-like structures are more stable than rings. As soon as the 
number of molecules is sufficient for the formation of a cage composed of fused rings whose 
(mean) molecular planes are either nearly orthogonal to each other (as in small (H20)„ 
clusters with n = 6 or 8) or form angles up to ca. 135° (as in (H20)„ clusters with n — 
12 or 20), such arrangement becomes preferable. A large number of cluster orbitals have 
mixed nature: they act as a bonding (or antibonding) orbitals in one plane (ring) and as 
TT bonding orbitals in the orthogonal subspacc. It is worth noting that in this orthogonal 
subspace, a continuous delocalized system, which spans several structural rings, is formed. 
The typical range of interplanar angles is determined by the mutual orientation of the initial 
molecular orbitals and their consistent contributions to o and tt binding. Moreover, even 
those orbitals, which look like antibonding in 0-H. . . O region, have such spatial character 
that act as bonding in 0-H and H. . . O subdomains. 

Thus, even in cage-like clusters, there are delocalized conjugated vr systems, which span 
either individual molecular rings or hoops composed of fused rings. (Recall that in the case 
of hydrogen fluoride, a sequence of quasi- aromatic segments rather than a continuous system 
can only be formed.) Besides that, tetracoordinated water molecules can be effective joints 
between molecular rings in a three-dimensional network. Hence, diverse structural groups 
can be present in water ensembles as long-living structure units. The strength and lifetime 
of an individual H-bond depend on its involvement in the conjugated tt system of a larger 
spatial fragment. Though it is a time consuming task, it is quite possible to determine 
all long-living structure fragments. However, it is already clear that fused molecular rings 
involved in closed or open hoops should have the largest lifetimes among long-living species 
due to general delocalized tt binding. 

4.5 Ionic clusters 

Character of the electron density redistribution that accompanies the formation of (HF)„, 
(H20)„, and (NHs)^ clusters is generally similar, and peculiarities are determined by the 
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degeneracy of molecular orbitals and the "part" of the electron density that is not involved 
in intramolecular binding. Then, it is interesting to consider ionic clusters, and primarily 
HF^ and HsO^, which are known to be stabilized by the strongest hydrogen bonds. 

In HF^, the binding is provided by electrons described by a a-orbital and two degenerate 
mutually orthogonal tt orbitals (Fig. [TTk shows one of them, yj^). Moreover, zeros of two 
corresponding antibonding tt orbitals accurately coincide with the proton position (Fig. [TTh 
shows one of them, iph)- Recall that in neutral molecular clusters, an alternation-sign surface 
could pass near a proton, but did not involve it. As a result, judging from the shape of the 
orbitals, one can say that there is some kind of a triple bond between two fluorine atoms. 
Moreover, both bonding and antibonding orbitals describe electrons that provide binding in 
F. . . H domains. 




Figure 11: Typical tt orbitals of (a) HFg and (b) HsO^. 



The quasi-triple F. . . F bond and the nearly zero electron density of antibonding orbitals 
in the vicinity of the bridge proton represent the limiting variant of the effects, which are less 
pronounced in neutral clusters. However, it is the effects that clearly explain the experimental 
NMR signatures of hydrogen bonds, namely, proton deshielding and spin-spin couphng of 
the nuclei of electronegative atoms. 

A similar picture (though more complex) is observed in the case of HsO^. In this ion, 
there are five cluster orbitals, each of which look like two embraced constant-sign segments 
(two of them are shown in Fig. [TTb . ipc and (fd)- Accordingly, each orbital describes electronic 
binding of two II2O fragments (with the zeros close to the position of the central proton). 
Thus, electrons of all five orbitals contribute to the binding in the system. As a result, the 
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bonds in H5O2 ion are energetically very close to those in HF^ ion. 

Being highly sensitive to the arrangement of nuclei, the peculiar electron density dis- 
tribution in these ionic structures survives only in the absence of external perturbation or 
when the perturbation has the same symmetry (no less than C2). Some examples of clusters, 
in which HF2 or H5O2 fragment steadily exists, are shown in Fig. [12] (there is much less 
freedom in constructing HF-based clusters). Note that in all systems, there is still a quasi- 
triple bond between two central oxygen or fluorine nuclei. However, for the spatial reasons, 
this kind of binding can by no means be spread on the residual structure. Therefore, all 
other hydrogen bonds are characterized by nearly the same electron density distribution as 
in neutral clusters. 



Figure 12: Examples of F (HF)„ and H+(H20)„ structures that involve stable HF2 and 
H5O2 fragments respectively. 

Thus, hydrogen atom plays the role of a bridge between two neighboring electronegative 
atoms, i.e., acts as a binding joint between their electron density "clouds". This role is more 
pronounced in the ionic systems, but could already be deduced from the electron density 
distributions in neutral clusters. If we consider a general situation, the formation of an 
X-H. . . Y-Z hydrogen bond is determined by the overlap of X-H and Y-Z bonding orbitals 
(or Y lone pair orbital); and hydrogen atom acts as a bridge. 





5 Conclusions 



The discovered peculiarities of the electron density distribution in the hydrogen bond do- 
mains enable us to analyze some criteria and characteristics of the bonds. 
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First of all, it is reasonable to note that there is always a and tt binding in H-bond 
domains. Hydrogen bonds are usually quasi-double and sometimes ( as in the case of HF^ 
and H5O2 ions) quasi-triple. 

H-bonds should be divided in two classes, namely, (i) separate bonds formed between 
two certain particles and (ii) bond networks that unite a large number of molecules and 
are stable when a delocalized conjugated tt system appears. Networks can be formed by 
hydrogen fluoride, water, and ammonia molecules, as well as any molecules that involve 
sterically accessible 0-H or N-H functional groups. The most stable building blocks of the 
networks are those with delocalized quasi- aromatic n-systems. 

Insofar as the proton in X-H. . . Y fragment acts as a bridge between the electron density 
"clouds" of X and Y, and the quasi-double or quasi-triple bond is formed between X and Y, it 
is natural that the bridge nuclei should he on the X. . . Y line. Then, the closer the X-H. . . Y 
angle to 180°, the stronger the bond. At the same time, the interaction is the strongest when 
X-H bond in the donor molecule and the orbital in the acceptor molecule that provides the 
main contribution to a binding have a common symmetry axis. This means that there is 
always a preferable orientation of the proton-acceptor molecule with respect to the H-bond 
direction, and the known directionality of hydrogen bonds becomes a natural consequence. 

The discovered peculiarities of the electron density distribution in H-bond domains clarify 
the origin of the NMR signatures of the bonds, particularly deshielding of H-bond protons 
and spin-spin interaction between X and Y. As we could see, electrons described by a and 
TT orbitals actually bind X and Y, especially when a conjugated system forms. At the same 
time, in many cases, zeros of both cr and tt functions are close to the position of the bridge 
proton and coincide with it in HFj and HsO^ ions. 

The existence of conjugated tt system in molecular clusters determines the cooperativity 
phenomenon. The spatial extension of cooperative effects depends on the extension of the 
delocalized conjugated n system, and the effects are most pronounced within the regions of 
quasi- aromatic ir-binding typical of structural rings and closed or open hoops composed of 
fused rings. 
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